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The dynamics of QCD is known to simplify somewhat in the limit where the number 
of colors, N, is large Q. Only planar graphs contribute to scattering amplitudes and the 
resulting perturbation series exhibits some of the qualitative features of the strong interac- 
tions. However so far, no explicit solution is available in the large N limit and it has thus 
led to very few quantitative results. 

Recently, Kazakov and Migdal || have proposed a novel approach which has the hope 
of providing an exact solution of the large N limit of QCD. They consider induced QCD 
which is obtained by integrating over the scalar fields in the lattice gauge theory with the 
partition function 
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where 4>{x) are NxN Hermitean matrices which reside on lattice sites x, U(xy) are unitary 
NxN matrices which reside on links < xy > between neighboring sites x and y, dip is the 
Euclidean integration measure for Hermitean matrices, [dU] is the invariant Haar measure for 
integration over the unitary group U(N) and V[<p] is a potential for the scalars. This model is 
invariant under the gauge transformations 4>(x) — > u(x)<fi(x)uj^(x), U(xy) — > u(x)U(xy)u :ii (y) 
where u(x) is an element of U(N). By restricting the trace of <fi to zero and the determinant 
of U to one in ([!]) we could also consider a model with SU(N) gauge symmetry. 

This model is soluble in the large N limit. By explicitly integrating over the unitary 
matrices in (|l|) (see Appendix A) one obtains a field theory for the eigenvalues of 4> 



J x i <xy> A[0(x)]A[0(y)] 



where A[0] = det^-^)- 7-1 = lli<j(0i — 4>j) is the Vandermonde determinant for 0. The 
eigenvalues 0; behave like a master field since the large N limit in @ is the classical limit 
and the integral can be performed by saddle point approximation. Migdal || has given 
an expression for the saddle point in terms of the density of eigenvalues. Corrections to 
the classical behavior and the spectrum of elementary excitations can also be computed [Q . 
This model has been considered further in U - [pjj] . 

If the model (f[J) has a second order phase transition and if the fluctuations in the vicinity 
of the critical point are non-Gaussian, one might expect that the critical behavior should 
be represented by QCD, the only known nontrivial four dimensional field theory with non- 
Abelian gauge symmetry. However it was pointed out in H that, like all adjoint lattice 
models, the Kazakov-Migdal model has an extra gauge symmetry which is not a symmetry 
of continuum QCD. The action in (|l|) is invariant under redefining any of the gauge matrices 
by an element of the center of the gauge group, U(xy) z(xy)U(xy), <p(x) — > 4>(x), where 
z(xy) G U(l) if the gauge group is U(N) and z(xy) G Z N if the gauge group is SU(N). 
(We shall call the symmetry a gauge symmetry in either case. This symmetry is the 
maximal subgroup of the transformations discussed by Gross [^0| and Boulatov [TJ| which 



can be implemented with field independent elements z(xy).) Because of this symmetry the 
conventional Wilson loop observables of lattice gauge theory have vanishing average unless 
they have either equal numbers of U and U' operators on each link or else, in the case of 
SU(N), unless they have an integer multiple of N IPs or iV U'*s. 
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In conventional QCD, the expectation value of the Wilson loop operator gives the free 
energy for a process which creates a heavy quark-antiquark pair, separates them for some 
time and lets them annihilate. From the asymptotics for large loops, one extracts the 
interaction potential for the quarks. If the expectation value of the Wilson loop behaves 
asymptotically like e~ aA where A is the area of a minimal surface whose boundary is the loop, 
the quark-antiquark potential grows linearly with separation at large distances and quarks 
are confined. The parameter a is the string tension. On the other hand if the expectation 
value of the Wilson loop goes like the exponential its perimeter then the potential is not 
confining. In the Kazakov-Migdal model, due to the symmetry, the expectation value of 
the Wilson Loop is identically zero for all loops with non-zero area. We can interpret this 
as giving an area law with a = oo, and no propagation of colored objects is allowed at all. 
(An exception is the baryon (U ) loops in the case of SU(N) where the correct statement 
is that A r -ality cannot propagate.) It is for this reason that the original Kazakov-Migdal 
model has difficulty describing pure gluo-dynamics. 

There are currently several points of view on how to cure this difficulty. In || it was 
suggested that if there is a phase transition so that the Zn symmetry is represented in a Higgs 
phase, the resulting large distance theory would resemble conventional QCD. This approach 
has been pursued in M, JTT]], |T2|. An alternative, which was advocated in 0,0], is to use 



unconventional observables such as filled Wilson loops which reduce to the usual Wilson 
loop in the naive continuum limit but which are invariant under Z^. The third possibility is 
to break the Zn symmetry explicitly. This was suggested by Migdal [|14j] in his mixed model 
in which he breaks the Zn symmetry by introducing into the model heavy quarks in the 
fundamental representation of the gauge group. In this Letter we shall examine a different 
explicit Ztv symmetry breaking and the use of filled Wilson loop observables. We shall show 
that the filled Wilson loops arise naturally from ordinary Wilson loops in a modified version 
of the Kazakov-Migdal model which has additional explicit symmetry breaking terms. We 
shall argue that one version of this modified model should be solvable in the large N limit. 
It thus provides a link between the approach which uses unconventional observables to solve 
the Zjy problem and the approach in which the symmetry is explicitly broken. 

We begin with a brief description of the filled Wilson loop operators which were intro- 
duced in || and discussed in detail in ||. These are a special class of correlation functions 
which survive the Zn symmetry. They are defined by considering an oriented closed curve 
T made of links of the lattice. The ordinary Wilson loop operator on T is given by 

Wqr] = Trj J] U{xy)\. (3) 
[<xy>er ) 

For any surface S which is made of plaquettes such that the boundary of S is the curve T 
we define 

w F [r,S\ = w[r] JJ w*[n] (4) 

□eS 

where □ denotes an elementary plaquette in the surface S. The filled Wilson loop for T is 
now defined as 

W F [T] = J2KS)W F [T,S] (5) 

s 
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where the sum is over all surfaces S whose boundary is the loop T with some (yet to be 
specified) weight function fi(S). Notice that for each plaquette □ G S we have inserted the 
negatively oriented Wilson loop VK'[n]. Thus for arbitrary weight functional n{S) the filled 
Wilson loop operator is invariant under the local gauge symmetry since it has equal 
numbers of U and operators on each link. Although we have assumed that the loop is 
filled with elementary plaquettes this can be easily generalized to other fillings (the other 
extreme case being the adjoint loop W[r]W^[r]). We can also define the 'filled correlator" 
of more than one loop by summing over all surfaces whose boundary is given by those loops. 

We now review some of the properties of the filled Wilson loop. In ref. |J it was shown 
that computing the expectation value of W0?[r] is equivalent to computing the partition func- 
tion of a certain statistical model on a random two-dimensional lattice. When computing 
Zn gauge invariant correlation functions of ?7-matrices in the master field approximation 
the ^-integral is evaluated by substituting the master field = diag(0 l5 . . . ,4>n) f° r the 
eigenvalues of 0. 

„ Jd4[dU]e*(Z*Wu iljl ...Uj lh --- 

If we consider for the moment surfaces which are not self-intersecting so that the filled 
Wilson loop correlator has at most one UW pair on any link we need to consider only the 
two field correlator < U^U^ >. Gauge invariance implies that 

f\dU]e N ^ Tr ( mu ^\U--\ 2 

< UuUL >= CMjk with Ca = ^ — l - ^ ,_ * J| (7) 

3 kl 3 3 3 J[dU]e N ^ Tl (' l>u ^ ) 

Thus, in the master field approximation, the expectation value of the filled Wilson loop is 
given by 



<w F [v] >=J2Ks) 

s 



( N \ 

11 Y, II C i{x) j{x 



i = itcs i(x)=l i !»«• 
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) 

links 



(8) 



This is a generalized Potts model on a random surface in which iV-component spins reside 
at each site and the Boltzmann weights Cy for the bonds are correctly normalized to be 
conditional probabilities; J2iCij=l, J2jCij=l. 

Techniques for evaluating Cij for general N and for arbitrary <fi are presented in [16 and 



jl7i . An explicit formula for SU{2) is given in ||. Although the general formula for in 
SU(N) is quite difficult to deal with, it is still possible to estimate the surface dependence 
of the statistical model partition function in (H) when <p is homogeneous by considering two 
different limits. (The details of the computation of CV, in these limits is described in the 
Appendix.) First is the limit in which is small. We call this the "High Temperature" limit 
since in this limit 
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is independent of i and j. It thus represents the Bolzman weights for a highly disordered 
system. We also consider the limit in which is large. We call this the "Low Temperature" 
limit since in this limit 

C£ = 5 lj + ... (10) 

and the value of the spin at each site is equal. In this case the represent the Bolzman 
weight for a perfectly ordered system. These two cases lead to profoundly different behavior 
for the filled Wilson loop. We shall assume that, by choosing the potential for the scalar 
field in (|2]) appropriately, either of these limits could be obtained (the eigenvalue repulsion 
due to the Vandermonde determinants in @ and the possibility of adding repulsive central 
potentials makes the low temperature limit more natural). 

We begin by estimating the value of the filled Wilson loop for a fixed surface S. In the 
"high temperature" case the statistical model is disordered. The sums over configurations 
at the various sites are independent and they contribute an overall factor N v (where V is 
the number of vertices on the surface) to the expectation value of the filled Wilson loop. 
Furthermore each link contributes a factor Cij=l/N so that the links contribute a total 
factor of N~ L where L is the total number of links. It follows that the expectation value of 
the filled Wilson loop goes like 

< W F [T, S] > HT ~ N y - L = N ^3(S) N -A(S) (n) 

where A(S) is the area and g(S) is the genus of the surface S (i.e. the number of plaquettes 
comprising S) and we have used Euler's theorem, x = 2 — 2g = V — L + A. We thus get the 
renormalization of the string tension 5anT = log N. Notice also that higher genus surfaces 
are suppressed and that the loop (genus) expansion parameter is 1/N 2 . This is precisely 
what is obtained in the conventional strong coupling expansion of Wilson's lattice gauge 
theory which is known to describe a string theory with extra degrees of freedom associated 
with self- intersections of the string p2|]. 

In the "low temperature" case, the statistical system is ordered. The spins on all the 
sites are frozen at a uniform value. In this case the partition function is proportional to the 
degeneracy of the ground state, 

< W F [T,S] > LT = iV (12) 

Note that in this case the statistical model gives no contribution to the string tension 
(<5a LT ~ 0) and there is no suppression of higher genus surfaces. 

In order to proceed to the evaluation of the filled Wilson loop we need to choose a weight 
function fi(S) in order to perform the sum over surfaces. The most reasonable criterion for 
choosing such a weight function is our desire to get a finite physical string tension in the 
continuum limit. In order to accomplish this goal we must choose a weight function n(S) 
which depend on the area of the surface differently in the low and in the high temperature 
cases. It is known that the number of closed surfaces with a given area grows exponentially 
as 

n(A) ~ A K{9) e^ A (13) 
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where n{g) is a universal constant which depends only on the genus of the surface and /xo is a 
non- universal, regulator dependent constant which will lead to a renormalization of the 
string tension. In our case, although the surfaces are open, the above formula should still be 
valid for surfaces whose area is much larger than the area of the minimal surface bounded 
by T. If the continuum limit of our theory is realized in the "high temperature" phase we 
should use the weight function (jl ht (S) ~ ]\[ A (s) e -^oA(S) _ This leads to a vanishing string 
tension in the lattice theory which is a necessary condition for having a finite string tension 
in the continuum limit. To accomplish the same goal in the "low temperature" phase we 
should use fi LT (S) ~ e -voA(S)_ Although these choices of fi(S) give the desired result, it is 
rather unnatural to have to choose fi(S) in such an ad hoc fashion. 

Fortunately there is a very natural way to obtain the sum over surfaces in (§). Consider 
the following expectation value 

where W[T] is the conventional Wilson loop. Remember that the average is weighted by the 
Kazakov-Migdal action as in (^): In the master field limit it is computed by integrating only 
over [/-matrices with = and with the Kazakov-Migdal action. Note that the exponent 
in ([TJ|) is simply the conventional Wilson kinetic term for the gauge fields in lattice gauge 
theory. If we expand the right hand side of (|14]) in A the non-vanishing terms are all of 
those surfaces which fill the Wilson loop. The result is thus a filled Wilson loop with a 
surface weight fi(S) = \ A ( S \ It is clear that we could obtain exactly the same expression (in 
the master field approximation) by evaluating the expectation value of the ordinary Wilson 
loop operator in a modified version of the Kazakov-Migdal model in which a conventional 
Wilson term (A^ D + is added to the action. This term breaks the Zn 

gauge symmetry explicitly and allows Wilson loop operators with non-zero area to have 
non-zero expectation values. We would expect that it is necessary to keep A small if one 
is to maintain the successes of the Kazakov-Migdal model. We shall now argue that in the 
"low temperature" limit this picture is self-consistent in the sense that the physical string 
tension is finite when A is small and consequently the saddle point solution of the original 
model is unchanged. We shall also see that this is not the case in the "high temperature" 
phase. 

Let us begin by determining how A should behave in the continuum limit if we are to 
have a finite physical string tension. As discussed above a necessary condition for having 
a finite physical string tension is that the string tension in lattice units should vanish. It 
is thus necessary for the bare string tension — In A to be chosen so as to precisely cancel 
the renormalization of the string tension due to both the statistical model to the sum over 
surfaces. It is straightforward to check that in the "high temperature" phase we must choose 
Aht = iVe~ Mo , whereas in the "low temperature" phase we must chooseA LT = e~ Mo . Notice 
that this A is proportional to N in the "high temperature" phase and thus cannot be assumed 
small in large N. 

In the large N limit of conventional lattice gauge theory the coefficient of the Wilson 
term must be proportional to N if one is to obtain a consistent large N expansion. In our 
case we see that this is true for the "high temperature" phase in which case the Wilson term 
is of the same order as the Kazakov-Migdal term and it thus plays an important role in the 
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infinite N limit. One can say that in this phase we have ordinary QCD. Unfortunately it is 
impossible to preserve the master field solution of the Kazakov-Migdal model in this limit 
since the Wilson term, being of order N, would modify the large N solution, ruining the 
self-consistency of the mean-field approximation as described here. 

The situation is much more appealing in the "low temperature" phase. In this case the 
required coefficient of the Wilson term is of order one. It is subdominant and therefore 
negligible in the large N limit. Thus, Migdal's solution |5| of the Kazakov-Migdal model 
in the large N limit should still apply to our proposed modification of the action. In fact 
the only reason that the Wilson term is important at all in the large N limit of the "low 
temperature" phase is related to the collective phenomenon which orders the statistical 
system on the surfaces. It effectively makes the statistical model's contribution to the string 
tension much smaller than would be expected from naive counting of powers of 1/N and a 
truly infinitesimal breaking of the gauge symmetry (A LT /iV — > as N — > oo) is sufficient 
to make the averages of Wilson loop operators non-vanishing. The self-consistency of this 
picture can also be demonstrated by computing the contribution of the Wilson term to the 
free energy. This can be computed in a small A expansion. For a cubic lattice the result is: 

Z=< e ME-^W^])) >= Zrm exp (W D(jP 2 ~ 1} (AL + 2A 6 LT + ...)) (15) 

is of order N (where V is the volume, D is the dimension). This should be compared with 
the free energy in the pure Kazakov-Migdal model which is proportional to N 2 . Here, the 
first term in the free energy is the contribution of the doubled elementary plaquette and 
the second term is due to the two orientations of the elementary cube. It is interesting 
that, to order 6, there is no energy of interaction of doubled elementary plaquettes with 
each other. We conjecture that the interaction energy of surfaces is absent to all orders and 
the free energy obtains contributions from all possible topologically distinct surfaces which 
can be built from elementary plaquettes. This suggests a free string picture of the "low 
temperature" limit of the Kazakov-Migdal model at lattice scales. 

We have thus far neglected the self-intersecting surfaces in the sum (B) which are gener- 



ated by the expansion of (14) in A. In order to evaluate the contribution of these surfaces 
we need to compute the correlator of n UU^ pairs on the same link. The computation of 
these correlators in full generality is quite complicated. In the Appendix we compute them 
in the "low temperature" (ordered) phase. We find that 

< U hjl . . . U injn Ul lh . . . U\ nln >= 5 hjl . . . 5 knln S£± (16) 

where S^ 1 '"^ is the tensor which is one if permutation of k± . . . k n and is zero 

otherwise. It is now evident that in this limit the ^/-matrices are replaced by unit matrices 
which freeze together the spin degrees of freedom on the various intersecting surfaces. As a 
special case we can consider a single, connected, self-intersecting surface. In this case all the 
spin indices on the surface are equal and since = 1 when all arguments are equal the 
partition function of the statistical model corresponding to that surface is simply N just as 
it was for a non-intersecting surface. Thus just as the statistical model does not contribute 
to the string tension it also does not contribute to the interaction energy of self-intersecting 
surfaces. This implies that in the "low temperature" limit, the sum over connected surfaces 
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which have a common boundary behaves like a Nambu-Goto string theory with no internal 
degrees of freedom. 

In summary, the self-consistency of the "low temperature" limit leads us to a new large 
N limit of the conventional lattice gauge theory coupled to scalars: 

d(f>[dU)expl-NY,TrV[(f>{x))+N ]T Tr(j)(x)U{xy)(f){y)U\xy) + 

\ x <x,y> 

+A^(D) + iy t (0))j (17) 

The conventional large N limit occurs when A is of order N and describes scalar QCD. 
The other limit occurs when N — > oo with A of order one. This model is soluble using the 
Kazakov-Migdal approach. 

It is the latter case in which A remains constant that is of special interest to us. In this 
case we saw that the large N expansion corresponds to a string theory with some unusual 
features. The partition function and the Wilson loop expectation value can be described as 
a sum over surfaces. What is unusual is that the genus of the surfaces is not suppressed in 
the large N limit, as it is in continuum QCD. (We do of course expect the higher corrections 
in 1/N to suppress higher genus terms.) For a continuum string theory this sum over the 
genus is badly divergent. This, together with the presence of tachyons, suggests that the true 
ground state of the string theory is some sort of condensate. This could pose a complication 
for the present version of the Kazakov-Migdal model in the continuum limit and deserves 
further attention. It is still a mystery to us how the sum over all surfaces at the lattice scale 
should turn into the sum over planar diagrams in the continuum theory of QCD. 

An alternative to the model presented here is the mixed model which was invented 



by Migdal [14| to solve the problem of symmetry. It contains heavy quarks in the 
fundamental representation of the gauge group. Despite the obvious differences between our 
model and Migdal's mixed model they have many features in common. As in all cases when 
there are fields in the fundamental representation, the asymptotics of the Wilson loops in 
the mixed model exhibit a perimeter law. In conventional QCD one would expect that if 
the quarks are heavy enough, there is an area law for small enough loops, i.e. there would 
exist a size scale which is far enough into the infrared region that the quark potential is 
linear but the interaction energy is not yet large enough that it is screened by producing 
quark- ant iquark pairs. Thus, in QCD we expect that adding heavy quarks would not ruin 
the area law for Wilson loops smaller then some scale. 

The mixed model has just the opposite scenario, it is possible to get an area law only when 
the heavy quarks are light enough. This is a result of the fact that, in the Kazakov-Migdal 
model, no Wilson loops are allowed at all unless the symmetry is explicitly broken. In 
the mixed model, the Z^ charge of links in a Wilson loop must be screened by the heavy 
quarks. This can happen in two ways. First, the Wilson loop can just bind a heavy quark 
to form an adjoint loop - giving a perimeter law for the free energy of the loop. This is 
the leading behavior if the fermion mass, M, is large. The free energy would go like 1/M P 
where P is the perimeter. The only way an area law might arise is when the fermions are 
light enough that their propagators could from a filled Wilson loop with free energy 1/M 2L 
where L « 2A is the number of links. Then, since the entropy for filled loops is much larger 
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than that for adjoint loops, these configurations would be important if M 4 < e M0 . Then, the 
asymptotics behavior of the Wilson loop would still have a perimeter law but there would 
be loops with 4A — P < /j, Q / InM where there would be an approximate area law. 
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APPENDIX: Asymptotics of Correlators 

The key to the solvability of the Kazakov-Migdal model is the fact that the single-link 
Itzyksen-Zuber integral can be done analytically p4| , |25 | 



[d U]e^^\ 2 = d f^^ (18) 



In this Appendix we consider two limiting cases of this integral and its correlators. These 



have already been solved in [ 16,17] . Here we present a quick and simple derivation in two 
limiting cases. 

1. Low Temperature: First consider the case where <fti are large. We also assume 

that the eigenvalues fa are not too close to each other in the sense that J2i^i 1 - ^ « N 

(4>i-4>j) 

(these two conditions are satisfied in the case of the semicircle distribution of the eigenvalues). 
The integral in (|I^) is known to be exact in the semi-classical approximation (see M for a 
discussion). The classical equation of motion is [C/0C/',0] = which, since <fi is diagonal, 
is solved by any U of the form Uq = DP where D is a diagonal unitary matrix and P is 
a matrix which permutes the eigenvalues, (P(j)P)ij = Sij(pp^y Also, when N is large and (f> 
is not too small the identity permutation gives the smallest contribution to the action in 
(|T8|) and therefore is the dominant classical solution. In this case we use this minimum to 
evaluate the correlators, 

Jfi / [dU\e N ^^U hh . . . U injn Ut lh . . . Ut nln = 6 hh . . . 5 knln S£X (19) 

where we have written the normalized integral over diagonal matrices 

S£X = lUid6 e U P<q sm 2 (6 p - ^)e l K + - +e -^^-^)//n£^n p<g sin 2 (^ - 9 q ) 

1 if i\ . . . i n is a permutation of k\ . . . k n , . 

otherwise 

We have decomposed the integration over unitary matrices into an integration over the 
diagonals and an integration over the unitary group modulo diagonals ||26|| . The diagonals are 
the 'zero modes' for the semiclassical integral and must be integrated exactly. The unitary 
modulo diagonal integral is damped by the integrand and is performed by substituting the 
classical configuration. Of course, to get the next to leading order the latter integration 
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must be done in a Gaussian approximation. It can be done for the first few correlators, The 
result is 

11 11 

Cu = 1 - — — —^2 when * ^ 3 Cij = 



k^i (0j - (f) k 

Similar calculations can be easily done for correlators of more than two C/'s 

i = j , k = £ Ca^k = 1 ^ ^ 
i^j , k=i 



1 1 

+ 



Oi - <Pn) 2 (<Pk - (/>n) 2 
1 1 



N (0, - 0,) 2 
11 1 



iv^ (0, - <^.) 2 (0 fe - <j, t y 

where the next corrections will be of the form 1 _ - 2 J2 m ^k 

\4>i-4>j) 

for the semicircle distribution of the eigenvalues one has 



(21) 



(22) 



Note also that 



1 



S (<i>j - <t>i) 



n — , for n > 2 



2. High Temperature: This corresponds to small 0. We obtain the correlators by Taylor 
expansion: 

C l3 = J[dU] (l + iVTr0f/0f/t + ...)| f/ij |2 = i_ + M + ... 



10 



REFERENCES 



[1] G. t'Hooft, Nucl. Phys. B72, 461 (1974). 

[2] V.A.Kazakov and A.A.Migdal, Induced QCD at large N, Paris / Princeton preprint 

LPTENS-92/15 / PUPT-1322, (June, 1992) 
[3] A.A.Migdal, Exact solution of induced lattice gauge theory at large N, Princeton preprint 

PUPT-1323 (June, 1992) 
[4] A.A.Migdal, 1/N expansion and particle spectrum in induced QCD, Princeton preprint 

PUPT-1332 (July, 1992) 
[5] I.I.Kogan, G.W.Semenoff and N.Weiss, Induced QCD and hidden local Zn symmetry, 

UBC preprint UBCTP-92-022 (June, 1992), Phys. Rev. Lett., in press. 
[6] I.I.Kogan, A.Morozov, G.W.Semenoff and N.Weiss, Area law and continuum limit in 

"induced QCD", UBC preprint UBCTP-92-026 (July, 1992), Nucl. Phys. B, in press. 
[7] A.Gocksch and Y.Shen, The phase diagram of the N = 2 Kazakov-Migdal model, BNL 

preprint (July, 1992) 

[8] S.B.Khokhlachev and Yu.M.Makeenko, The problem of large-N phase transition in 
Kazakov-Migdal model of induced QCD , ITEP-YM-5-92, (July, 1992) 

[9] M.Caselle, A. D. 'Adda and S.Panzeri, Exact solution of D=l Kazakov-Migdal induced 
gauge theory, Turin preprint DFTT 38/92 (July, 1992) 
[10] D.J. Gross, Some remarks about induced QCD, Princeton preprint PUPT-1335 (August, 
1992) 

[11] Yu. M. Makeenko, Large-N reduction, Master field and loop equations in Kazakov-Migdal 

model, ITEP-YM-6-92 (August, 1992) 
[12] S. Khokhlachev and Yu. Makeenko, Adjoint Fermions Induce QCD, ITEP-YM-7-92 

(August, 1992). 

[13] I.I. Kogan, A.Morozov, G.W. Semenoff and N. Weiss, Continuum limits of "induced 
QCD": lessons of the Gaussian model at D = 1 and beyond, UBCTP 92-27, ITEP- 
M7/92 (August, 1992), Int. J. Mod. Phys., in press. 

[14] A.A.Migdal, Mixed model of induced QCD, Paris preprint LPTENS-92/23, (August, 
1992). 

[15] A. Migdal, Phase Transitions in Induced QCD, LPTENS-92/22 (August, 1992). 

[16] S.L. Shatashvili, Correlation functions in the Itzykson-Zuber model, IAS preprint 

IASSNS-HEP-92/61, (September, 1992). 
[17] A. Yu. Morozov, Correlators in the Itzyksen-Zuber integral, ITEP-M10/92, (September, 

1992). 

[18] D. V. Boulatov, Local Symmetry in the Kazakov-Migdal Model, NBI-HE-92-62 (Septem- 
ber, 1992). 

[19] B. Russakov, From Hermitean matrix model to lattice gauge theory, TAUP 1996-92 
(September, 1992). 

[20] A. Migdal, Bose condensation and Z N symmetry breaking in the mixed model of induced 

QCD, PUPT-1343 (September, 1992). 
[21] D. V. Boulatov, Infinite-tension strings at d>l, NBI-HE-92-78 (November, 1992). 
[22] V. A. Kazakov, Phys. Lett. 128B, 316 (1983); Zh. Eksp. Teor. Fiz. 85, 1887 (1983) [Sov. 

Phys. JETP 58 1096 (1983). 
[23] A. B. Zamolodchikov, Phys. Lett. 117B, 87 (1982). 
[24] C.Itzykson and J.B. Zuber, J. Math. Phys. 21 (1980), 411. 



11 



[25] Harish-Chandra, Amer. J. Math. 79 (1957) 87. 

[26] Despite the simple form of S its tensor structure is quite complicated. For example 

^feifea _ °k 1 °k 2 + °k 2 °k 1 ~ °i2°k 2 



12 



